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Abstract 

Effects of electron correlation on the Fermi surface is investigated for the two- 
dimensional Hubbard model by the quantum Monte Carlo method. At first, an 
infinitesimal doping from the half filling is focused on and the momentum dependent 
charge susceptibility nik) = —j^- is calculated at a finite temperature. At the 

temperature T ~ jj, it shows peak structure at (zb7r/2,zb7r/2) on the Fermi surface 
(line) . It is consistent with the mean- field prediction of the <i-wave pairing state or 
the staggerd flux state. This momentum dependent structure disappears at the high 
temperature T ~ U. After summarizing the results of the half filling case, we also 
discuss the effects of the doping on the momentum dependent charge susceptibility. 
The anisotropic structure at half filling fades out with sufficient doping. 
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Effects of electron-electron interaction have been a subject of intense investi- 
gation in condensed matter physics. One of the most remarkable phenomena 
is the Mott transition, which is a quantum phase transition driven by the 
interaction. At a rational filling, the ground state of the interacting electron 
systems sometimes belongs to the Mott insulator. Furthermore, various kinds 
of instabilities toward, e.g. superconductivity or antiferromagnetism exist near 
the Mott transition. The electronic properties near the Mott transition can 
deviate from a simple Fermi liquid and it leads to interacting phenomena. 
In particular, the behavior of the quasiparticles near the Fermi surface is an 
important topic to be investigated. 

Recently singular momentum dependence of the low energy excitations has 
been focused on in the strongly correlated electron system. In this context, 
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the spectral weight of the two-dimensional Hubbard model has been investi- 
gated [1-3]. Moreover, the deformation of the Fermi surface due to the interac- 
tion is observed in the t-J model [4,5]. Also the angle resolved photoemission 
spectroscopy (ARPES) experiments suggest anisotropic properties in the low 
energy excitations [6]. 

We investigate the momentum dependent charge susceptibility n(k) = 
in the two-dimensional Hubbard model by the Quantum Monte Carlo (QMC) 
simulations. The n(k) = dn j^ measures a charge fluctuation with momentum 
resolution. Since the infinitesimal shift of the chemical potential \x corresponds 
to an infinitesimal doping at a finite temperature, the n(k) also shows the dis- 
tribution of an infinitesimally doped carriers in momentum space, 5n(k) ~ 
K,(k)5fi. Without interaction, the n(k) is peaked uniformly on the Fermi sur- 
face (line). At half filling, it a square in the Brillouin zone, \kx\ ~\~ \ky\ — tt. 
Even in the presence of the interaction, we call this square the Fermi surface 
(FS) at half filling in this paper. 

The ground state of the two-dimensional Hubbard model at half filling belongs 
to an insulating fix point (antiferromagnetic Mott insulator). Since the Mott 
insulator has a finite charge gap, the charge fluctuation does not exist at zero 
temperature. However, it is interesting to discuss the effects of the interaction 
on the charge fluctuation at a finite temperature. We study correlation effects 
on the momentum dependent charge susceptibility the n(k), focusing on the 
difference between (n,0) and (tt/2, 7r/2) points. Many kinds of fixed points can 
play a role in the charge properties at a finite temperature even if they can 
not be the long-range ordered at the zero temperature. Therefore, comparing 
the QMC results with possible mean-field solutions is interesting [7]. After 
summarizing the discussion at half filling, we also study the effects of finite 
doping on the re(fc) in this paper. 

The two-dimensional Hubbard model on a square lattice is 
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where (j, k) denotes nearest-neighbor links, U the Coulomb interaction, and 
t hopping amplitude. Periodic boundary condition is imposed. We obtain 
approximation-free results by using the finite temperature auxiliary field quan- 
tum Monte Carlo (QMC) method [8-10]. To calculate the momentum depen- 
dent charge susceptibility K(k), we directly evaluated the n{k) in the QMC 
simulations as follows, 
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where (3 denotes an inverse temperature and L the linear system size. 

At half filling, the sign problem does not occur due to the particle-hole sym- 
metry. Since the simulation is carried out in the grand canonical ensemble, 
information on the infinitesimally doped systems are statistically taken into 
account, even though the averaged system is half-filled. Away from half fill- 
ing, we focus on the relatively large doping case to avoid the difficulty of the 
negative sign problem. 

At first, let us discuss the half filling case [7]. Figure 1 shows the n(k) for 
U/t = and U/t — 4 at T/t — 0.2. Without interaction, the n(k) is peaked 
on the FS and its value is constant there. On the other hand, with the inter- 
action U/t — 4, the re(fe) shows peak structure around (±7r/2,±7r/2) points. It 
gives an anisotropic structure on the FS and indicates that the (±7r/2,±7r/2) 
points are more sensitive to the shift of the chemical potential than (±7r,0) or 
(0, ±7r) points. One may say that the holes are doped near the (±7r/2,±7r/2) 
points upon a small doping. Figure 2 shows the temperature dependence of 
the n(k) for U/t — 4 at several temperatures. The peak structure at the 
(±7r/2, ±7r/2) points is clear at T ~ t 2 /U and becomes ambiguous as the tem- 
perature increases and finally vanishes above T ~ U . Here we note that, at 
these temperatures, the antiferromagnetic correlation length is smaller than 
the linear system size and the size dependence is negligible. Our results suggest 
that the interaction brings about the anisotropy in the charge fluctuation. The 
anisotropy is evident at the temperatures of order of the effective exchange 
interaction J ~ t 2 /U. For the high temperatures T > U where the Coulomb 
interaction U becomes irrelevant, the anisotropic feature vanishes [7]. 

Although the ground state belongs to the antiferromagnetic Mott insulator, 
a different picture can appear as a crossover at an intermediate temperature 
T ~ J. Comparing the QMC results with possible three different mean-field 
solutions ( rf-wave pairing state, staggered flux state, and Neel state), the 
peak structure is qualitatively reproduced by the ci-wave pairing state or the 
staggered flux state [7]. 

Next, we discuss doping effects on n(k). Due to a severe negative sign problem, 
it is difficult to treat systems with very small doping (near the half-filled) . The 
result of the n(k) for the doped system is shown in Figure 3 with the filling 
(J2k n (k)) /L 2 ~ 0.64. The correlation effects suppress the charge fluctuation 
and the sum of the n(k) over the Brillouin zone is substantially reduced. 
Further the n(k) is spread over the total Brillouin zone for the interacting case 
(U/t = 4) in comparison with the non-interacting case (U/t = 0). However, as 
for the anisotropy on the Fermi surface (line), a clear peak structure does not 
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observed. This is clearly different from the case of the half filling. It implies 
that, with sufficient doping, instabilities which is driven by the interaction may 
become irrelevant and the deviation from a simple Fermi liquid is negligible. 
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(a) U/t=0 




Fig. 1. The momentum-resolved compressibility «(fc) = dn(k)/dfi obtained by the 
quantum Monte Carlo simulations at T/t = 0.2 on a 16 x 16 lattice: (a) for U/t = 0; 
(b) for U/t = 4 [7]. Without interaction ([//i = 0), the n(k) is constant on the 
Fermi surface. On the other hand, for U/t = 4, the re(fe) shows peak structure at 
(±7r/2, ±7r/2) points on the Fermi surface. 




Fig. 2. The momentum-resolved compressibility re(fe) = dn(k)/dn obtained by the 
quantum Monte Carlo simulations on a 16 x 16 lattice for U/t = 4 at various tem- 
peratures [7]. The peak structure at (±7r/2, ±7r/2) points vanishes at temperatures 
T > U. 
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(a) U/t=0 




Fig. 3. The momentum-resolved compressibility fc(fc) = dn{k)/dji obtained by the 
quantum Monte Carlo simulations at T/t = 0.2 on a 12 x 12 lattice: (a) for U/t = 0; 
(b) for U/t = 4. The filling is set to be 0.64. In contrast to the half filling, a clear 
peak structure does not observed. 
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